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ABSTRACT 



The objective of this investigation is to derive the reliability and the 
associated probabilistic failure criterion for composite materials under 
combined stress. In the analytical derivation, the concept of joint probability 
was used and applied to the Weibull distribution function. In applications, given 
the experimental measurements of the necessary statistical parameters for the 
specific composite, the probabilistic criterion of the composite failure and the 
reliability of the specific structure can be predicted. 

Graphical representations for the joint reliability and joint failure contours 
were made in two and three dimensional space for the several different sets of 
statistical strength parameters to illustrate the effect of parameters on reliability. 
Such understanding will enhance selection of fiber and matrix (which have their 
own statistical strength parameters) and can lead to improvement in reliability of 
some composite components in an aircraft structure. These reliability and 
failure concepts can also be used in repair problems by selecting the proper 
composite material with the appropriate statistical parameters. 
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I. INTRODUCTION 



A. FAILURE AND RELIABILITY OF THE COMPOSITE MATERIAL 

In a composite, the matrix and reinforcement materials are combined to form 
a macroscopically multiphase material . The macroscopic state gives rise to a 
material which is internally redundant against local defects and damages. 
[Ref.1] The ply is the basic unit in constructing composites which consists of 
fibers embedded in a matrix. It is possible, therefore, to relate composite design 
properties to corresponding constituent properties. As a result, it is possible to 
design structural components made of fiber-composite when the design 
requirements and the properties of the candidate constituent materials are 
known. For high performance structures, the unidirectional fiber-matrix 
composites consist of strong stiff fibers embedded in a comparatively ductile 
matrix . For a rational structural design, the strength of such unidirectional ply 
needs to be expressed in terms of the stress tensor. Because the ply is 
composed of numerous fibers, the composite failure process is sequential, 
therefore ad hoc failure criteria can not be applied to analyze the composite 
failure. Furthermore, such conventional failure criteria such as Mises and 
Tresca's criteria are related only to a single parameter (mean function of 
failure). As a result, they can not be applied to multiple parameter functions, so 
suitable general relations to represent failure contours of composite materials 
are required. 

In the application of a one-dimensional failure criterion, the magnitude of 
stress (one component) at each spatial location of the structure is mapped into 
respective points on the stress space (which is a line). If all the points within the 
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domain are bounded by the tensile strength limit and the compressive strength 
limit, then the entire structure may be considered safe. The reliability under 
unidirectional stress, then, is the product of the reliability of each spatial 
segment at the respective stress level. Under planar combined stress, the 
magnitude of stresses (3 components) at each spatial location of the structure is 
mapped into respective points on the six dimensional stress hyperspace. So the 
stress that is interior to a domain bounded by the failure surface is safe and the 
failure criterion is the analytical, graphical, or numerical representation of the 
failure surface. To formulate the failure criterion, consistency of mathematical 
operations can be assured by adhering to the quotient rule of tensor variables. 
Historically, the formulation of many failure criteria were intuitively taken to be 
based on failure mechanisms but in fact, they are phenomenological because 
of the inconsistency between macro and microscopic failure modes. 
Phenomenological failure criterion may be considered as a mathematical 
model of the material transfer function relating the external excitation (stresses) 
to the material's response (failure). The phenomenological mathematical 
model is intended to aid experimental design to facilitate interpolation, 
correlation, and retrieval of experimental observations. They do not, in general, 
address the statistical variability of strength. 

For a combined stress failure mechanism, when one stress component does 
not affect the strength of another component, the failure is mechanistically 
uncoupled: conversely when one stress component affects the strength of 
another component, the failure is mechanistically coupled. In addition, if 
probabilistic failure is also taken into consideration, the condition for the failure 
mechanism is not fixed; it depends on the stochastic combinations of the 
intrinsic strengths. Under combined stress, mechanistic uncoupling with 
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statistic independence can interact to produce phenomenologically coupled 
contours. On the other hand, mechanistically coupling with statistical 
dependency can also produce phenomenologically coupled statistical contours. 
These cases can only be distinguished by examination of both the contour in 
the stress space and the statistical strength distribution in the probability space. 

In any phenomenological failure criterion, the shape of the failure envelope 
is never completely known until experiments are performed for all possible 
combined states of stress. As a result, complete experimental determination in 
two or more dimensions place an impractical demand on time and resources 
because of the large number of tests needed. However, with the establishment 
of a probabilistic failure criterion, this experimental failure contour can be 
related to the shape parameters (the variability) in the material’s principal 
direction. With such supplemental information, failure contours connecting the 
mean value can be established by a smaller number of samples. Therefore an 
analytical model can complement the failure characterization by reducing the 
required number of experimental measurements, providing bases for 
comparison of the material’s performance, and simplifying data and data base 
reduction. In structural analysis applications, the reliability of a specific structure 
under a combined stress can be determined by mapping the spatial domain into 
the stress domain through stress analysis. Using the probability failure criterion, 
the reliability of each spatial location corresponding to the respective stress 
tensor and magnitude can be calculated. The joint reliability of the individual 
spatial location can in turn be combined to predict the reliability of the entire 
structure. 
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B. OBJECTIVE OVERVIEW 

The objective of this investigation is to derive an analytical model for 
composites failure and reliability under combined stress conditions. The 
Weibull distribution function was used in the probabilistic modeling. Chapter 2 
presents the theoretical formulation for failure and reliability, and the resulting 
equations are shown graphically in the following chapters. Chapter 3 shows 
the joint reliability and linearized failure Cumulative Distribution Function (CDF) 
in two and three dimensional space. The effects of parameters which depend 
on the material were analyzed. Chapter 4 shows the joint failure CDF in two 
and three dimensional space. The effects of material parameters were also 
analyzed. In Appendix B, these probability distributions were compared to 
experimental data availabe in the literature. 

Finally, Chapter 5 provides conclusions and some remarks for the future 
work in this area. 
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II. THEORETICAL FORMULATION 



A. RANDOM VARIABLES 

A failure criterion for a solid prescribes the risk of failure in terms of the 
causes of the failure. Mechanical catastrophic rupture is considered herein. 
The effect of mechanical failure may be considered to be caused by either 
stress or stain depending on the microscopic failure processes. If the 
constitutive relation of the material is available, conversion between stress and 
strain can be readily made and the two variables may be considered as 
equivalent. Similarly, a probabilistic failure criterion which describes the 
probable risk of failure can be expressed in terms of stress or strain as random 
variables. Stress is selected as the random variable herein with the 
understanding that the variable can be changed to strain with the appropriate 
constitutive relation. We denote the probabilistic random variable as Xj, for the 
scalar component of the stress tensor , and Xj for the realized random variable 
(i.e. the strength Xj under the respective stress component q). 

B. PROBABILISTIC CONDITIONS OF FAILURE 

Under combined stress conditions, when more than one component of the 
stress tensor assumes non-zero values, the effect of each component of the 
stress tensor contributes to the failure of the composite. Failure of the 
composite does not occur if failure by each individual stress component does 
not occur. 

For the composite subjected to a planar state of stress, we denote the failure 
events as: 



5 



A; Failure caused by normal stress along the fiber direction - a-] 

B: Failure caused by normal stress perpendicular to the fiber direction - 02 
C: Failure caused by shear stress along the fiber direction - 03. 

Then the failure of the composites is either by A or by B or by C or by all three as 
illustrated in Figure 2-1: 

P(A + B -h C) = P(A) + P(B) + P(C) - {P(AB) -h P(BC) -t- P(CA)} + P(ABC) 

( 2 - 1 ) 

As a special case, if all three events are mutually exclusive, then Equation 
(2-1) can be simplified by 

P(A + B + C) = P(A) -h P(B) -t- P(C) (2-2) 




Figure 2-1 : Joint Probability 
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C. JOINT DISTRIBUTION FOR FAILURE 



1 . Discrete Case 

In the development of the analytical expressions for the probability of 
failure, we first developed the general form of the distribution for three random 
variables using the notation Xj for the random variables and xj for the realized 
random variables. For the planar stress, the customary notation for the range of 
i = 1 , 2,6 is used. These general expressions are then specialized for the 
combined stress cases by replacing the realized random variables by Oj. 

The probability of failure under the combined influence of three random 
variables can be expressed by the joint probability density function (pdf) of the 
individual realized random variables: 

P(X, = o,, Xj = = Og) = f (a,, 



where f(a-| , 02, 05) > 0 and 'LYL f(a-| , 02 > ~ fandom variable X0 is to 

take on any one of the values of 051 , oq2, <753,..., 05^, , then the joint pdf at 
some combination of Xi=a-|j, X2=a2j, and X0=G0|^ can be expanded by 

P(X, = o,i,Xj =Oj.,X, = Oe^) = 



n 

P(X, = 0,.,X2 = Oj ) =f(0,i,0j ) = If(0,i.02,.06k) 

' k=0 

I 

= °2i'^6 = "6k) = *("2r "6k) = 

^ i=0 



( 2 - 5 ) 



( 2 - 6 ) 



j=0 



( 2 - 7 ) 



and the marginal probability function can be obtained for the three random 
variable case: 
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( 2 - 8 ) 



m n 



P(^ = <’li) = = 2; If (o,^.a 

j=Ok=0 ^ ^ 



I n 



P(^2=°2i)= '(®2i) = iSIf(0,i.a2j,0e,) 



i=0k=0 
I m 



P(^6 = <’6k) = '(<’6k) = ^2;f(o,|,a ..a^A 

i=0j=0 ^ ' 



(2-9) 



( 2 - 10 ) 



which can be written 

I rn 



If,(o,i) = 1. If 2(0 ) = 1. If = 1 

j=0 ^ ^ k=0 



i=0 



I m n 

IIIf(<5,2<J2,OgA=1 

k0j=0k=0 ^ ' 



(2-11) 



( 2 - 12 ) 



2. Continuous Case 

If three random variables are considered, then Equation (2-4) through 
(2-7) can be expanded by 

P(X, i c,, Xj £ Cj, £ Ce) = F(a,, Oj, Oe) 



a a a 
6 r 2 r 1 



= J I j f (u,v,w)dudvdw 



— 00 — 00 — 00 



a o 

2 r 1 



P(X^<o^,X2 < 02) = F^a^.a^) = j J jf (u,v,w)dwdudv 



— 00 — 00 — 00 



o a 
6 2 



P(X2 <a2,Xg < Cg) = F(^a2,ag) = J J Jf (u,v,w)dudvdw 



— 00 — 00 — 00 



^6 f 



P(X^<a^,Xg < Gg) = F(a^,CTg) = J J Jf (u,v,w)dvdudw 



— 00 — 00 — o 



(2-13) 



(2-14) 



(2-15) 



(2-16) 



and the marginal distribution can be obtained by 

O 

,| 00 00 

P(X^ < = F(a^^ = J J Jf(u,v,w)dvdwdu 



— 00 —00 —00 



(2-17) 
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( 2 - 18 ) 



P(X2 < Gg) = F^O^) = J I Jf(u,v,w)dwdudv 



— oo — oo — oo 



P(Xg < Gg) = F(Og) = I I Jf(u,v,w)dudvdw 



a 

^ 6 



( 2 - 19 ) 



which can be written 




G2 given that c-\ occurs can be represented by a function form: 




( 2 - 21 ) 



where f(G-| , G2) = P(X-| =g-| , X 2 =G 2 ) 
a. General Case 

If we consider the combined stress case, the random variables X-| , 
X2, and X0 can be assumed to be the failure strength in the 1 , 2, and 6 direction 
respectively. And under a combined stress condition, the failure occurs when 
either failure strength X-| , X2, or X0 is reached. Therefore the probability of 
failure under combined stress X-| and X2 can be represented by 



and then the above relation can be related to the joint distribution function as 
follows 



F(gi , G2) = P(Xi < Gi , X2< G2) = P(Xi + X 2 : Xi < g^ , X 2 < G 2 ) (2-22) 

and similarly for three random variables case: 

F(g-|,G2,G0) = P(X-| <G-| ,X2^G2,X0<G0) =P(X-|+X2 + X0 : 



P(Xi + X2) = P(Xi ) + P(X2) - P(Xi , X2) 



X-| <G-| ,X2^G2>X0<G0) 



( 2 - 23 ) 
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Equation (2-22) is proven in Appendix A-1 using the Weibull 
distribution function for the independent case. 

If the joint distribution of failure under the combined stresses of a-\ , a2> 
and a0 is considered, the strength failure criterion of a composite is the 
minimum of the failure strength X-| , X2, and X0, that is, at X^ = X-|aX2aX 0. ( We 
are not dealing with the minimum magnitude among X-|, X2, and X0 but rather 
one value among X-| , X2, and X0 which is less than the intrinsic failure 
strengthes of , 02, and oq becomes the minimum value here.) So 
F(a^,a2,ag) = P{Xc = X^ aX^ a 

where ( a-| , a2< ct 0 ) is the three dimensional stress domain. If the relations 
among a-| , 02 and gq are given by 

ai =Bi2a2, ^2 = B26(76> = (2*24) 

then, the joint pdf P(X-|< a-\ ,X2^G2 ,Xq<gq) can be evaluated. If a0 is assumed 
to be zero, then the joint probability functions, P(X-|<a-| ,X2^02) can be divided 
into two probability distribution functions corresponding to S-| and S2 domains 
as shown in Figure 2-2. The domain Si is applied for the case when 
[Xi / X2 ^ B-|2 and X-| < a-|] which the composite fails by failure strenth 
X2 = X-|AX2, whereas the domain S2 is applied for the case when 
[X1/X2 ^ B-|2 and X2 < 02] which the composite fails by failure strength 
X-| = X-|AX2. So the following joint relation can be obtained: 

P(X-| < Gi , X2 ^ G2) = P(X-| /X2 ^ B-| 2 and X-| ^ g-| ) 

+ P(Xi /X2 > Bi 2 and X2 < 02) (2-25) 

Applying Equation (2-25) to domains described in Figure 2-2, then 
P(X^<o^,X2<a2) = J Jf (u,v)dudv-i- J Jf(u,v)dudv 

Si Sg 
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Using end limits and considering P(Xi<a-| ,X 2 ^ci 2 ) = F(a-| ,a 2 ),then the joint 
distribution function for the general case can be represented by 



^2 <• 

F(ai,a2) = J J_2_f(u,v)dvdu+ ^f(u,v)dudv 

- ''12 



0 B 



12 



(2-26) 



where u and v are dummy variables. 




iL 

02 



Figure 2-2 : Domains For The Joint Distribution 



b. Specific Case 

The joint probability of the bi-axial case (Equation (2-25)) can be 
expanded for specific cases. 

Independent Case . If the events X-j = a-| and X 2 = <^2 
independent for all a-| and © 2 , then the joint probability can be represented by 
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(2-27) 



P(Xi=ai,X2=a2) = P(Xi=a-|) P{X2=a2) 
or equivalently 

f(ai,a2) = fl(a-|) f2(o2) 

If three independent random variables are considered: 
P(Xi=ai.X2=a2,X6=06) = P(Xi=ai) P(X2=02) P(X6=06) 
f(a-| ,02,00) = f(oi ) f(o2) f(o0) (2-28) 

Similarly, if the random variable Xi and X 2 are independent for all 



o-| and 02 ' 

P(Xi<oi,X 2<02) = P(Xi<Oi) P(X2<02) 
or equivalently; 

F(oi,02) = Fi(oi)F2(o2) (2-29) 

where Fi(o-|) and F 2 (o 2 ) are the marginal distribution functions of X-| and X 2 
respectively. Similarly, for the three independent random variable case: 

F(Xi<oi,X 2<O2,X6<O0) = P(Xi<oi) P(X2<02) P(X6<Oe) 
or equivalently 



F(oi ,02,00) = F(oi ) F(o2) F(o0) (2-30) 

So the joint distribution function can be expanded a step further for 
the independent case. 

From Equation (2-26) (refer to Appendix A-2) 




U >1 



du-J \(v)Fi(B^2'')^v 



(2-31 a) 



Equation (2-31 a) can be further expanded by integration by parts 
= ^(^l) + ^" 2 (^ 2 ) - 2Fi(C^0^2(^2) 



12 o V 127 o 



(2-31 b) 
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Comparing Equations (2-31 a) and (2-31 b) reveals that the 
arguments of the pdf and CDF are exchanged so either Equation (2-31 a) or 
(2-31 b) can be used depending on simplicity of algebra. 

(a) Weibull Distribution Model. If the Weibull function is applied to 
the independent case, then Equation (2-31 a) can be expanded by 



F(<J,,02) = 



a 



1 - exp] - 






pj JJ 



1 - exp' 



V^2 y 



a 

2 



-J 



a -1 

V U A ^ 



P,APi 



a 



_U 

pj j 



1 - exp] - 



“ 2 I 






du 



cc 1 

-J ' 



a 



2 y 



^P2 



exp^ - 






1 - exp 






P, 



du 



By rearranging terms, the joint distribution function for the 
Weibull model can be obtained by 



’ An 



a. 



a. 



U 



'2rct 



^ If f 



a “1 
2 



exp 



^ 12^2 
a 



du 



a 



U 

IPs 



P, 



fdu 



(2-32) 



( 2 ) Independent And Identical Case . Equation (2-31 a) can be further 
expanded if the identical case is considered. 

For the identical case, F(a) = F-| (a^ ) = ^ 2 (^ 2 ^ equivalently f(o) 
= fi(G-|) =f 2 (o 2 ), so expanding Equation (2-31 a) 

= ''(''1) + '"(" 2 ) - fo MF(B,2V)dv 



and similarly Equation (2-31 b) can be represented by 



(2-33a) 
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(2-33b) 



F(a-i .G2) = F(ai ) + F(a2) - 2 F(oi ) F(a2) 

For the special case when B 12 = 1 . <? = o-] =02 Equations 
(2-33a) and (2-33b) can be simplified to 

F(o,,a2) = 2F(o) - 2jf’ (u)F(u)du 



and rearranging finally gives the joint distribution function for independent and 
identical case. (refer to Appendix A-3): 

F(oi,02) = 2F(o)-(F(o))2 (2-34) 

Equations (2-33a) and (2-33b) can be evaluated for specific 
models such as the Weibull distribution function. 

(a) Weibull Distribution Model. For the Weibull model, the pdf (f(x)) 
and cumulative density function (CDF: F(x)) can be represented by 



"•’'■(tKr-'Hs)'! 

pi -(«)") 



Substituting Equation (2-35) and (2-36) into Equation (2-33a) 



F(o,,o^) = 



“I 



1 expj 



1 -exm - 



^2^ 



ct'l 



J JJ 



C'(Hir"|‘{-(sK 


1 - exp 


1 (v)l 


1. (fXf) “'{-(Sll 


1 - exp“ - 
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By integrating and rearranging the terms, the joint distribution function of the 
Weibull model under bi-axial stress a-| and 02 can be simplified to (refer to 
Appendix A-4) 



^(^1’^2) = I 






a 









(2-37) 



D. MEAN FOR THE JOINT DISTRIBUTIONS 



1. General Case 

If X-| and X 2 are two continuous random variables having the joint 
probability density function f(a^,a 2 ), the mean or expectation of X-|and X 2 can 
be obtained by: 



— 00—00 
00 00 

^ —00 —00 ' 



(2-38a) 

(2-39a) 



Considering the stress, the bottom limits of the double integration are 
changed, so Equations (2-38a) and (2-39a) can be expanded specifically for 
the stress case which we are concerned with 



00 00 



(2-38b) 

(2-39b) 



Similarly if three random variables are considered, then Equations 
(2-38b) and (2-39b) can be expanded by 

00 00 00 

>^1 = Jo Jo Jo'’l'{°v'^2’°6)'‘<’l'''"2'‘<’6 



00 00 00 



t*2 = Jo Jo Jo°2'(°1'°2'°6)'‘"i'‘''2'''^( 



(2-40a) 

(2-40b) 
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(2-40C) 



0 0 0 ^ ■' 



2. Specific Case 

a. Independent Case 

If the random variables X-| and X 2 are independent for all a-| and 02 

f(ai,G2) = fi(a-|)f2(a2) 

and Equations (2-38b) and (2-39b) can be expanded by 



00 00 

^‘2 = i„ ij,'’2*l('"l)'2(‘’2)'''’l'''’2 



(2-41 a) 
(2-42a) 



taking the integral and noting that F-|(oo) = F 2 (~) = 1, then Equations (2-41 a) 
and (2-42a) can be simplified to 



^‘1 = 



J^o,f,(o,)da, 



(2-41 b) 



^^2 = 



J„°2*2('’2)''°2 



(2-42b) 



Similarly, for the three random variable case: 

00 00 00 00 

" -^0 -^0 -^0 



(2-43a) 



00 00 00 



-^0 '’2'l('’l)*2('’2)'6('’6)'^°,'‘‘’2^°6=i„''2*2(°2)'^°; 



00 00 00 



(2-43b) 



‘*6 = 4 '’ 6 ',('',)* 2 ('’ 2 )' 6 ('’ 6 )'‘°,'''" 2^°6 = Ns ' 



(2-43C) 

Observing Equations (2-41 b) through (2-43c), it can be concluded 
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that the mean values of the random variables X-| , X 2 , and X 0 are dependent 
only on their corresponding probability density functions. 



b. Independent And Identical Case 

For the independent and identical case, the joint probability function 
can be obtained by 

f(oi ,02,05) = f(ai)f (02) f(06) 
where f(o) = f-| (o) = f 2 (o) = f6(o) 

Applying these relations to Equations (2-41 b) through (2-43c), the 
mean values for the joint distributions can be simplified for the independent and 
identical case. 



Ml Weibull Distribution Model . If the Weibull probability function is 
considered. Equation (2-44) can be expanded as follows 



0 

Using the definition of the Gamma function, we finally get the 
mean function: 



E. JOINT RELIABILITY 
1. Concept Of Reliability 

The common notion of reliability is the confidence in the ability of a 




(2-44) 




letting = (o/p)^, then 





(2-45) 
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device to perform adequately over a range of specified conditions. Such 
qualitative measurement of confidence and adequacy can be quantitatively 
calculated by probability. In the formalism of probability, the conditions over 
which the device is intended to perform is defined as the random variable X: 
and the probability over the range of conditions is defined as the CDF over the 
range that the random variable experiences in service. That is, over the sample 
space a experienced in service, the probability of failure is defined by the CDF 
over the sample space: 

P(X < a) = F(a) 

Since reliability is the compliment of failure, it is defined by 

R(o) = 1 - F{a) (2-46) 

2. Joint Reliability 

If X-| and X 2 are two random variable, then the joint reliability of Xi and 
X 2 can be defined by 

R(oi ,02) = 1 - P(Xi <oi ,X2<02) (2-47) 

or equivalently for the three random variable case 
R(oi ,02,06) = 1 - P(Xi<oi,X2<02,X6<06) 

For the discrete case, equation (2-47) can be expanded by: 

~ S S(u,v) 

(2-48) 

FV 0^,02, Og) = 1 - X X X^(u.v,w) 

u<o.v<o.,w<o . 

1 2 6 



If the continuous case is considered, the joint reliability function of X-| and 
X2 is defined by 



o o 

F^0i,02) = 1 -J 7 ^Uu.v)dvdu 



(2-49) 



and from Equation (2-49) the marginal reliability function can be obtained: 
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(2-50a) 



R^(a^) = 1 -j Jf(u,v)dvdu 

— oo — oo 

^ 2 (^ 2 ) l^(u.v)dudv 



— 00 — 00 



(2-50b) 

If the three random variable cases are considered, Equation (2-50a) and 
(2-50b) can be further expanded: 

00 00 

= 1 - J j Jf (u,v,w)dvdwdu 



— 00 — 00 —00 



00 00 



*^2(‘^2) “ ^ “ i* I jf (u,v,w)dwdudv 

— 00 — 00 —00 

a 

0 00 00 

= 1 - j J jf (u,v,w)dudvdw 



— 00 — 00 — 00 



(2-51 a) 
(2-51 b) 
(2-51 c) 



Rl(oi), R2(<72)> R6(06) equivalent to the one dimensional 

reliability function which are defined by 



= ■' - J 



\(u)du 

» 

a 

^f 2(u)du 

- 00 
^6 

f g(u)du 



(2-52a) 

(2-52b) 

(2-52C) 



a. General Case 

The general equations for reliability can now be applied to the 
calculation of the mechanical strength reliability of composites under combined 
stresses by specifying the random variables as Xj and the realized random 
variables as cj. The joint reliability under combined stress a-| and 02 can be 
obtained by 
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where Xi and X2 are the failure strength in each coordinate(see Figure 2-2). 

Noting that the joint pdf under combined stress a-| and 02 is 
represented by the double integration of f(o-| , 02) over the domains S-| and S2 
and that the double integration of f(aj, 02 ) over the domain S3 (see Figure 2-2) 
is defined as the joint reliability under combined stress and 02 , then the 
following relation can be obtained: 



00 00 

R(a,,02) = 1 - F(0,.02) = f (u.v)dudv 



(2-53) 



Comparing Equation (2-53) to Equation (2-26), Equation (2-53) can 
be represented in another form: 

°1 00 °2 “ 

FVoo) = 1 -J J Lj f(u,v)dvdu-i- J Jf(u,v)dudv 

0 B~ 0 g 

12 ^12'^ (2-54) 

If the joint reliability under three combined stress is considered, 

equation (2-53) can be expanded by 

00 00 00 

= = Jo j^f(u,v,w)dudvdw 



6 2 1 



(2-55) 



b. Specific Case 

The joint reliability of bi-axial stress case can be further expanded for 
the specific case as in the probability distribution function. 

f1 f Independent Case . If the random variable X-| and X2 are 
independent for all o-j and 02 , then Equation (2-53) can be expanded: 



f\<'v''2) = Jo l„f,Mf2(V>dudv 

2 1 



(2-56) 



After integration, Equation (2-56) can be simplified to (refer to 
Appendix A-5) 

R(ai,a2) = {1-Fi(ai)}{1-F2(a2)} (2-57) 

noting that Ri (<?•)) = 1-F-|(a-|) and R2(<^2) = 
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R(ai ,02) = Ri (oi ) R2 (cj2) (2-58) 

For the three random variable case, Equations (2-56), (2-57), and 
(2-58) can be obtained by 

oo oo oo 
6 2 1 



= (1-Fi)(1-F2)(1-F6) 

= Ri (<71 ) R2(<72) Fieioe) (2-59) 

In application to a composite, if the probability of failure model (Fj) 
is known, then the reliability of the composite can be calculated. The probability 
of failure model can be inferred from the physical consideration of the failure 
processes. 



When a filament composite is loaded along the fiber direction, 
local failure begins when the weakest fiber fails. The load carried by the broken 
fiber is transferred to the neighboring fibers. Upon additional increase in load, 
additional fibers fail, leading to the increase of failure sites distributed over the 
composite. The higher the load, the higher the density of such failure sites and 
the higher the probability of clustering. The spatial clustering of the fiber failure 
sites leads to stress concentration and ultimately causes the catastrophic failure 
of the composite. Harlow and Phoenix investigated the probabilistic modeling 
of the above sequential failure events and arrived at a modified weakest link 
model in which the link of the chain is a bundle. [Ref. 2] For a limited range of 
the random variable, it can be approximated by the two parameter Weibull 
model [Ref. 3]: 






R(a,) = exp|-^^J I 

When a filament composite is loaded perpendicular to the fiber 
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direction, the failure mode is similar to that described by classical fracture 
mechanics. The exception is that for a composite, all inherent flaws are not 
randomly oriented, but are aligned with and propagate along the fiber. In such 
a case, the largest crack dominates and it is effectively a weakest link process 
over the physical volume. On these physical grounds we may use the Weibull 
model for both the transverse strength and shear strength: 




' 


f”’' 


a 
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I _ 




a ^ 
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> 




Ip,. 
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^^2y 





(2-60) 



Equation (2-57) can be specified for these Weibull distribution 

functions: 

equivalently for the three random variable case: 

FVa ,a ,a ) =exp^ 

^ ^ (2-61) 
If we note F(a-| , 02 ) = 1 - R(cji , a 2 ), then the joint distribution 
function obtained in Equation (2-32) can be simplified by 
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a ^ 
6 


< _ 




(^ 2 ) 
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IpJ 








F(a^,a 2 ,ag) = 1 -expj- 
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0 
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0 
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f-^1 
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IpJ 
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(2-62) 



(2-63) 



(2) Independent And Identical Case . If the random variables X-| and 
X 2 are independent for all a-| and 02 and the probability distribution function for 
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each random variable is the same, then F(a) = F-| (a-j ) = F 2 (o 2 ). Similarly R(o) = 
R-| (ai) = R2(<^2) ^he Equation (2-58) can be simplified to 

R(ai.a2) = R(ai)R(a2) (2-64) 

As a specific case, if the Weibull distribution function is 
considered. Equation (2-64) can be further expanded: 



Fi(o^,a 2 ) =exf 



f 




a 
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f ^2^ 
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r 


u. 


J - 


.T, 


I 



(2-65) 



For B-|2 = 1. then a = a-j = 02 and the joint reliability function for 
Weibull model is obtained. 



R(o,o) = exp{-2(^|)“} 



( 2 - 66 ) 



For the three random variable case. Equations (2-64) through 
(2-66) can be expanded by 



R(ai ,02,00) = R(oi) R(o2) R(oe) 

{ a 

-( p, 

R(o,o,o) = exp|-3^-|j 



a 



\ P J 



(2-67) 

( 2 - 68 ) 

(2-69) 
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III. GEOMETRIC REPRESENTATION OF JOINT RELIABILITY IJNDFR 

COMBINED STRESS 



A. BACKGROUND FOR GEOMETRIC REPRESENTATION 

CDF and pdf are useful functions which contain all the relevant information 
about the statistical properties of a random variable, in our case, the strength of 
composite. These functions are required in reliability analysis, design, 
acceptance, maintenance, and operational logistics. 

In applications, in order to map out the statistical failure surface even for the 
bivariate case (the biaxial combined stress), an exceedingly large number of 
experiments are required. A large number of experiments are frequently 
impractical due to economic and time constrains and other considerations. 
Therefore, it is important to visualize the shape of the failure surface in order to 
narrow the range of experiments to be focused on the critical regions in the 
stress domain. 

In this chapter we present the geometric representations for the reliability 
functions for the independent case derived in Chapter II. The graphical 
presentation is based on the probability plots and the failure surface 
representations, which were also investigated. Examination of these graphical 
representations will shed light on the appropriate experiments necessary for the 
identification of whether the failure processes are independent or dependent. 
When the independence is established, the entire failure surface (all 
permutations and combinations of combined stress case) can be calculated 
from the uniaxial strength statistics. 
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R ( s i g m a ) 



B. JOINT RELIABILITY FUNCTION 
1. Two Random Variable Case 



From Equations (2-56) and (2-66), the reliability function for the Weibull 



model can be obtained by 



R(a) = exp 




(3-1) 



RELIABILITY 




□ alpha=3 
♦ alpha=25 



Figure 3-1 : Effects Of a In Reliability 
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If we consider a/(3 as a normalized variable, then the reliability function 
versus a/(3 can be plotted for various values of a. Figure 3-1 shows two different 
graphes for large and small values of a. For the small values of a, the graph 
shows a relatively smooth curve, but as a increases, the graph approaches to 
the step function and these relations can be further expanded for the joint 
reliability problem. If the independent combined stress cases are considered, 
the joint reliability of the composite under combined stress can be obtained from 
Equation (2-58). 

R(ai,a2) = Ri(cJi) R2(cJ2) 

Expanding R(a-i , 02 ) 
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where 



R(Oj) = exp 
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|,R( 02 ) =exp' 






(2-60) 



If we assume that 

Pi = V 12 P2 

and noting that 

= Bi2 c^2 

then the joint reliability for the independent case can be represented by 

R(ai,02) = Ri(ai) R2(o2) 



(3-2) 



= exp 



a 












exp 






(3-3) 



If we further assume C-| 2 =^ 1 2 /^ 12 - ^^en Equation (3-3) can be simplified 



by 



26 



F^a^.a^) =exp‘ 

From Equation (3-4), the reliability functions Ri(oi) and R2 (ct 2) can be 
plotted independently in two dimensional space provided that a-\, a 2 , and C12 
are given. Then Ri(a-|) and R2(c?2) are combined to yield the joint reliability 
function of R(a-|,a2) as shown Figure 3-2. Observing the Figure 3-2, it can be 
noted that the joint reliability is affected by the smaller value between R-|(ai) 
and R2 (ct 2) a given value of ai/Pi , that is, the value of the joint reliability is 
always close to the smaller value between R-i(a-i) and R2 (ct 2)- When the values 
of both R-|(ai) and R2(r^2) ^re close to ’1’ for the specific range of a-i/Pi, the 
value of the joint reliability is also close to 1, but as either Ri(gi) or ^ 2 (^ 2 ) 
decreases, the joint reliability also decrease depending on the smaller value 
between R-| and R2. We can also observe the effect of C12 for the given a-| and 

« 2 - 

As shown in Figures 3-2a and 3-2c, the R2(o2) shifts to the right from the 
R-|(gi) as C12 decreases from ’1’ whereas R2(o2) shifts to the left as C-| 2 
increases from T’ while R-| (g-| ) remains constant. So, the joint reliability follows 
the smaller function between R-| (g-| ) and R2 (ct 2) sod appears to almost overlap 
when C-|2 is far from '1'. As the value of C-12 approaches to T, the values of 
R-i(g-i) and R2(<?2) approach each other and when C-|2=‘>. Rl(cfl) sod R2 (ct 2) 
cross each other at Gi/Pi=1. The joint reliability function exists to the left of 
Ri(g-|) and R2(o2) and has a weaker reliability than R-|(gi) or R2(o2)- So the 
joint reliability is always less than other uni-axial directional reliabilities. 




« n 



exp 



- C 



12 






(3-4) 
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R(sigma) 



RELIABILITY (A1=5,A2=3,C12=0.5) 




sigma1/beta1 



Figure 3-2a ; Reliability Vs (Ci 2=0-5) 



Q Ri 
♦ R2 
+ R 
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R ( s i g m a ) 



RELIABILITY (A1=5,A2=3,C12=1) 




sigma1/beta1 



Figure 3-2b : Reliability Vs Oj/Pi (Ci2=1) 



Q R1 
♦ R2 
+ R 
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R(sigma) 



RELIABILITY (A1=5,A2=3,C12=5) 




□ R1 
♦ R2 
+ R 



sigma1/beta1 



Figure 3-2c : Reliability Vs (Ci2=5) 
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Figure 3-3 shows the joint reliability in three dimensional space for 
different values of C-| 2 and for fixed a-] and 0.2- 




Figure 3-3 ; Reliability In 3-D Space (aj=5, a2=3) 
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Here we introduce the concept of linearized failure CDF to show the 
effect of the joint failure CDF at the tail area. 

To get the linearized failure CDF, denoted by F*(a), use logarithmic 
algebra. If we take the logarithm of Ri(a-|) and note that Ri(a-|) = 1-Fi(ai), 
then 



a 



In 



{R,(o,)} = ln{1-F,(o,)}=- 



vP,. 



(3-5) 



Multiplying both side by (-1) and taking the logarithm again, the 
linearized failure CDF can then be defined as 



(3-6) 



Substituting Equation (3-5) into Equation (3-6), then 



F,*(o,) = a.ln 






P. 



V 1 y 



(3-7) 



By the same way, linearized failure CDF corresponding to 02 t>e 
obtained: 



F/(o^) = a. 



ln(C,J + ln 






pjj 



(3-8) 



Comparing Equation (3-8) to (3-7), it can be noted that F 2 *(o 2 ) is 
vertically shifted and rotated when compared to F-|*(ai). This result is 
represented in Figure 3-4. 

Similarly F*(ai . 02 ) can be obtained by 



F*(a,,a^) = ln j 



fc 



+ 



c — 
S2P 



(3-9) 



Comparing Equation (3-9) to (3-7) and (3-8), it can be noted that 
F*(a-|,a 2 ) is no longer a linear relation in terms of ln(ai/Pi). Figures 3-4a, 3-4b, 
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and 3-4c show the linearized failure CDF for different values of Ci2- As shown 
in these graphs, F 2 *(cj 2 ) is shifted upward and rotated clockwise as C-| 2 
increases and this causes a higher probability of failure. The joint linearized 
failure CDF, F*(ai,a 2 ), then follows the higher value between Fi*(ai) and 
F2 *(cj 2) and does not show the linear line, especially near the area where 
F-|*(a-i) and F 2 *(a 2 ) intersect each other. So we can observe that , as the value 
of C-|2 increases for the specified a-| and a 2 , F2*(o2) shifts upward and rotates 
clockwise causing the intersection between Fi*(ai) and F 2 *(a 2 ) to go 
downward, or in other words creating a higher probability of failure. 



LINEARIZED FAILURE CDF 
(A1=5,A2=3,C12=0.5) 




Figure 3-4a : Linearized Failure CDF (C-i 2=0-5) 
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(sigma) 



LL 



LINEARIZED FAILURE CDF 
(A1=5,A2=3,C12=1) 




Q F1 
♦ F2 
+ F* 



In(sigma1/beta1) 



Figure 3-4b : Linearized Failure CDF (Ci2=1) 
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(sigma) 



LINEARIZED FAILURE CDF 
(A1=5,A2=3,C12=5) 




In(sigma1/bela1) 



Figure 3-4c : Linearized Failure CDF (C-|2=5) 
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R(sigma) 



As a special case, if a=ai=a 2 , then the joint reliability function and 
Equations (3-7) through (3-9) can be simplified by 






ra V 






0-0 






.P.y 



F,_*(a,) = a{ln(CJ + lr|^J. 

F*(Oj,o \ =alnO| + !n(l + c “,3 
V 1 J 



(3-10) 

(3-11) 

(3-12) 

(3-13) 



RELIABILITY (A1=A2=5, Cl 2=0.5) 




sigmal/betal 



Figure 3-5a : Reliability Vs ai/Pi(aj=a2=5) 



a R1 
♦ R2 
+ R 
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R(sigma) 



RELIABILITY (A1=A2=15,C12=0.5) 




□ R1 
♦ R2 
+ R 



sigma1/beta1 



Figure 3-5b : Reliability Vs CT|/Pj (a2=a2=15) 
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Comparing Equations (3-12) and (3-13) to Equation (3-10), it can be 
noted that F2*(o2) is shifted vertically by a ln(v-| 2/612)- F*(ai,02) is also shifted 
vertically by ln{1 -t-(vi 2/B1 2)} regardless of values of Ci 2- Figures 3-5 though 
3-7 show the reliability in two and three dimensional space and linearized 
failure CDF in two dimensional space for a = ai = a2- 




Figure 3-6a : Joint Reliability In 3-D (apa2=5) 
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Figure 3-6b : Joint Reliability In 3-D (ai=a2=15) 
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As shown in Figure 3-7, each failure CDF is parallel to each other and 
the joint linearized failure CDF almost coincides with F-|*(o-i) or F2*(a2), 
depending on the value of C 12 , and is parallel to both F-|*(a-|) and F2*{o2). So 
in the physical sense, C ^2 very important for estimating the failure and 
reliability. The joint reliability is affected by the weakest reliability function and 
this is true for the joint linearized failure CDF. 



LINEARIZED FAILURE CDF 
(A1=A2=5,C12=0.5) 




□ Ft* 
♦ F2* 

-I- F* 



In(sigma1/beta1) 



Figure 3-7a : Linearized Failure CDF {ai=a2=5) 



40 



LINEARIZED FAILURE CDF 
(A1=A2=15,C12=0.5) 




In(sigma1/beta1) 



Figure 3-7b: Linearized Failure CDF (apa2=15) 



R (sigma) 



2. Three Random Variable Case 



If the three random variable case is considered, the joint reliability under 
the combined stress can be obtained from the Equation (2-79): 

R(ai , <52, c5q) = Ri (ai ) R 2 (a 2 ) Re W 
Expanding this equation for the Weibull model: 
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( 2 - 61 ) 



RELIABILITY (A1=25,A2=5,A6=3) 
(C12=0.5,C16=1) 




Q R1 
♦ R2 
n R3 
+ R 



sigma1/beta1 



Figure 3-8a : Reliability Vs (Ci 2=0-5, C-|e=1) 
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R(sigma) 



RELIABILITY (A1=25,A2=5,A6=3) 
(C12=1,C16=1) 




0 12 3 
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Figure 3-8b : Reliability Vs (C-j 2 =^ 1 6=1) 
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If we assume that 



Pi = Vi2 p2. P 2 = ^16 P6 



and noting that 



(3-14) 



Oi = Bi 2 02, Oi = B10 05 



where 



B B = V V =— ^ 

12 B ’16 B ’12 V ’16 V 
21 61 21 61 



(3-15) 



(3-16) 



then the joint reliability for the independent case can be represented by 




Figure 3-9 : Joint Reliability In 3-D (Ci 0 = 1 ) 
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R(oi,a2,a6)=Ri (oi)R2(o2)R6(o6) 
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(3-17) 



where 



C, = 



V 



1 2 



12 B 



•C,e = 



V 



1 6 



B. 



(3-18) 



12 "16 

Furthermore, if C-|2 C 10 are given, then C 26 can be calculated 

using the chain rule: 



C„ = 



V 



26 






16 



26 B, 



'26 ^ 2^6 ^ 2^1 Pi Pe ^12^16 ^12 (3-19) 

So, if a-|, tt 2 , uq, C-| 2 . acd C -|6 are given, then the joint reliability 
function, R(a-| , 02 , 03 ), can be plotted in two dimensional space with respect to 
01 /Pi- Here we can also introduce the concept of linearized failure CDF to show 
the effect of the failure CDF at the tail area. To get the linearized failure CDF, 
denoted by F*(o), use the same procedure described in Equations (3-5) through 
(3-9), then 



'^2*(°2)=“2{'"(C,2)+'"[i;]} 

'"6'(°6)=“6{"XC,6)+'"[^]| 



(3-7) 

(3-8) 

(3-20) 
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(3-21) 



So it can be noted that the joint linearized failure CDF, F*(a-| .02,05) does 
not show a linear relation in terms of ln(o-|/Pi) except when 01=02=05, instead 
it is affected by Fi*(oi), F2*(o2), and F5*(o5). That is to say, F*(oi ,02,05) is 
dominated by the weakest value among Fi*(ai), F2*(o2), and F5*(o5) so 
F*(oi ,02,05) is always shown on the top of each of graphs as shown in Figures 
3-1 Oa and 3-1 Ob. Figures 3-8a through 3-8b show the reliability for different 
values of C12 when 01, 02, C15 is fixed. These figures are expanded to three 
dimensional space as shown in Figure 3-9. 



LINEARIZED FAILURE CDF (C12=0.5,C16=1) 
(A1=25,A2=5,A6=3) 




□ F1* 

♦ F2* 

n F3* 
+ F* 



Figure 3-1 Oa : Linearized Failure CDF (Ci 2=0.5, Ci5=1) 
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(sigma) 



LINEARIZED FAILURE CDF (C12=C16=1) 
(A1=25,A2=5,A6=3) 




In(sigma1/beta1) 



o F1* 
♦ F2* 

□ F3* 

+ F* 



Figure 3-1 Ob : Linearized Failure CDF (Ci 2 =^ 1 0=1) 
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As a special case, when a=a-|=a 2 =a 0 , then the joint reliability and 
linearized failure CDF functions can be simplified by 



Observing Equations (3-22) through (3-24), it can be noted that F-|*(ai), 
F2*(<^2)' F*(ai,O 2 ,cT 0 ) are parallel each other. And the Joint 

linearized failure CDF, F*(ai ,G2>cf6)' 's also affected by the largest function 
among Fi*, F 2 *, and F 0 *. The parameter a is related to the rotation whereas 
Cl 2 srid C 10 are related to the vertical shift of the graph. 




(3-22) 




(3-23) 
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IV. PROBABILISTIC FAILURE CONTOUR UNDER COMBINED STRESS 



A. PROBABILISTIC FAILURE 

For a combined stress failure mechanism, when one stress component does 
not affect the strength of another component, failure is mechanistically 
uncoupled, as shown in Figure 4-1. When one stress component affects the 
strength of another component, failure is then mechanistically coupled, as 
shown in Figure 4-2. 




Figure 4-1 : Mechanistically Uncoupled Mechanism 
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Figure 4-2 : Mechanistically Coupled Mechanism 



If the probablistic failure is considered, the mechanistically uncoupled case 
may contribute to statistical coupling so both failure mechanisms are not 
interacting and the probablistic failure mechanism is distinguished from the 
deterministic case. In order words, the condition for the failure mechanism is 
not fixed, it depends on the stochastic combinations of the intrinsic strengths. 
So, for the combined stress case, uncoupled independent statistical effects 
produce phenomenologically coupled statistical contours and the mechanically 
coupled dependent statistical effects also produce phenomenologically coupled 
statistical contours. In the following subsection, the joint probabilistic failure 
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CDF will be discussed for the two and three random variable cases and the 
effects of the parameters will be shown graphically. 



B. Two Random Variable Case 

To show the joint failure probability, we here introduce the Weibull 
distribution function. From Equation (2-62), the joint failure distribution function 
can be obtained for the two random variable case. 



From the above equation, (a-|/Pi) and { 02 /^ 2 ) computed if the values 

of F(c-j,a 2 ) and C 12 are given. Figure 4-3 shows the joint failure CDF in three 
dimensional space and Figure 4-4 shows corresponding joint failure contour in 
two dimensional space for the specific value of F. As shown in Figures 4-3 and 
4-4, the shapes of the failure function and failure contour depend on the 
parameter a, which is material dependent. For a small value of a, the failure 
contour shows a smooth curve but as the value of a increases, the joint failure 
contour approaches to the shape of a rectangle. This phenomenon is the same 
for the joint reliability contour as shown in Figure 4-5. It can also be noted that 
the failure contour is affected by a large value of a between a-| and a 2 - So if 
anyone of the values of a is larger when compared to another value of a, then 
the failure contour approaches to the shape of a rectangle, depending on the 
magnitude of a. These relations are shown in Figures 4-6 and 4-7. Another 
result, shown in these failure contour graphs, is that the distance between the 
contour line is small when a is large compared to those when a is small. So the 
intersection of the function on both axis depends on the value of each a. 




(2-62) 
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Joint Failure CDF 
0.00 0.29 0.58 0.87 




Figure 4-3a : Joint Failure CDF In 3-D (a^=a2=5) 
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Joint Failure CDF 
0.00 0.29 0.58 0.87 
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A1=A2=25 




Figure 4-3b : Joint Failure CDF In 3-D (apa2=25) 
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Figure 4-4a : Joint Failure Contour (ai=a2=5) 
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Figure 4-4b : Joint Failure Contour (ai=a2=25) 
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Figure 4-5a : Joint Reliability Contour (ai=a2=5) 
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Figure 4-5b : Joint Reliability Contour (aj=a2=25) 
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The contour line of F(ai,a 2 ) =0.5 exists around 1.0 in both the normalized x 
and y axis regardless of a, but as a increases, the F=0.5 contour line 
approaches to 1.0 in both normalized axis. 




Figure 4-6a : Joint Failure CDF (ai=5,a2=3) 
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Joint Failure CDF 

0.29 0.58 0.87 




Figure 4-6b : Joint Failure CDF (ap25,a2=3) 
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As an application, if we compute the stress of a specific part of system such 
as a rudder or ailerons using the finite element method, then we can estimate 
the probabilistic failure and reliability of that specific part of the system. 

The parameters a and p, which were used in the failure CDF and reliability 
function, are material dependent so these values should be determined through 
experimentation and then these parameters can be applied for that specific 
material. 




Figure 4-7a : Joint Failure Contour (ai=5,a2=3) 
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Figure 4-7b : Joint Failure Contour (ap25,a2=3) 
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Joint Failure CDF 



C. THREE RANDOM VARIABLE CASE 



If the three random variable case is considered, the joint failure CDF can be 



developed using Equation (2-63). 



a 



a 



a 






F(a a ,a ) = 1 -expj- - - 

^ I IPJ IP2J 










Figure 4-8a : Joint Failure CDF (C-|0=1) 
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Figure 4-8b ; Joint Failure CDF (Ci6=3) 
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From the above equation, normalized axis coordinates can be computed if 
the values of F(a-| . 02. oq), Ci 2, and Ci 5 sre given. And if we specify the value 
of C-10, then we can plot the joint failure CDF in three dimensional space with 
respect to a-|/Pi and 02/^2 shown in Figure 4-8. As the value of C-| q 
increases, the possibility of failure increases, that is to say, the intersection on 
the normalized x-axis moves to the left because 05/^5 term is embedded in the 
a-|/p-| axis. Figure 4-9 shows the corresponding failure contour line in two 
dimensional space for the specific value of F(a-| , 02, 05). 




Figure 4-9a : Joint Failure Contour (Ci0=1) 
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As proven graphically in the three random variable case, the probability of 
failure increases as the number of random variables increase. So in a physical 
sense, a specific composite material will show the highest reliability when the 
external loads are applied in the pricipal axis direction because there are only 
two random variable, X-\ and X 2 , in the reliability fuction. But when the external 
loads are applied in the off-principal axis direction, there exist some shear force 
random variable, X 0 , which decreases the reliability of the material. 




Fig. 4-9b : Joint Failure Contour (Ci2=3) 
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As an application, if the stress of the maneuvering aircraft is computed in an 
arbitrary direction, then the stress can be transformed into fiber, shear, and 
transverse matrix directions which can then be used to predict the probability of 
failure at that specific load. Sample calculations were made and the results 
were analyzed with respect to experimental data in Appendix B. 
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V. CONCLUSIONS 



This study was directed towards deriving the joint reliability and joint failure 
CDF under a combined stress condition. The results were specialized for the 
Weibull distribution function based on the observations that the composite 
material failure is adequately represented by the weakest link model. The effect 
of the statistical strength parameters, a and p, which are material dependent are 
illustrated using two and three dimensional graphical representations. 

To analyze the reliability of the composite, the inherent statistical strength 
parameters in the composite's principal direction need to be experimentally 
measured. Substitution of these parameters in the probabilistic failure criterion 
will allow for the estimation of the reliability of the composite for any state of 
stress. 

Comparison of joint failure distribution under the restriction of independence 
to experimental data suggests that mechanistic coupling of the failure 
mechanism needs to be included in future extensions of the formulation and 
data of much larger number of samples has to be performed at critical stress 
ratios to conclusively examine probabilistic independence. 

Further studies may be extended to apply these reliability and failure CDF to 
the specific part of a system which is made of a composite material. To do this, 
the stress of the specific part of a system such as an elevator or ailerons should 
be analyzed using such as finite element method. These equations can also be 
used in a repair problem which requires the least cost and most effective 
method to analyze the combined stress in a specific part of a system. 
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So, if we analyze the combined stress, the appropriate fiber and matrix with 
proper parameters can be selected to make the composite fit to that specific 
part. 
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APPENDIX A 



Chapter II Equations 



1 . For the independent case: 

P(X-|+X2 : X-|<ai, X2^G2) = P(Xi<a-|) + P(X2^a2) - P(X-|<ai) P(X2<a2) 

= Fi (ai ) + F2(a2) - F-| (ai ) F2(c2) 

If the Weibull distribution function is considered, Fi(a-|) and F 2 (a 2 ) can be 
substituted into the Weibull distribution function: 



P(X^ + X2:X^<a^,X2<a2) = 



1 - exp" 






an r 
1 



. I ^1 J J . 



1 - exp' 



a 



v^2y 



1 - exp' 






1 



i IP, 

By expanding and simplifying 



1 - exp' 



.P2 



a ^ 
2 



P(X^ + X^:X^ < g^,X2 ^^ 2 ) = ^- 

which is the same result as the Equation (2-32) 
2. For the independent case: 

a 







a 

1 




a 

2 




f OA 




f 0 ^') 




< 


1 




2 










<^2> 





I 12 J 

Taking the integrals inside the parenthesis first and substituting limits upon 
evaluation: 
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F,(o) = j'’f,(u)du,F2(<j) = /'f2(u)du 



Noting F-| (oo) = P 2 ('«) = 1 from the definition of CDF and expanding 



a a 

= l Mu)du-J f,(u)Figa-\ju 






+ J 'f,(v)dv-J 'f,(v)F/B v)dv 

o o 

Taking integrals and noting that F-j (0) = F2(0) = 0 

''(°V'’2) = ^(°l) + '^2(<’2) - 4 

a 

-I fj(v)F,(B v)dv 

o 

3. From Equation (2-33a) 

F(o a -) = 2 F(o) - 2 jt (u)F(u)du 

By integration by parts 



^(^ 1 , 02 ) =2F(a)-2 



Jf(u)F(u)du 



noting 

J^(u)F(u)du = 0.5{F(a)}^ 

O 

then the joint distribution function for independent case can be obtained by 
F(oi , 02 ) = 2 F(o) ■ {F(o))2 
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4. By expanding and integrating 

F(0,,02 




1 -expi - 



^2^ 






P 



r r“ ^ 
^12 

B“ + 1 



o 

^ 1 




a 



exfX-[-^l - 






- exps - 




1 ^ 



exp 



o 

2 






“ ''5,2'' 



P 



By substituting and expanding 






1 - expi - 1 



ct") n r 



+ 



expl- 



expi- 






P 






- 1 



n V 



1-expi-l-y^ 

P 






exp^ 



y ) j 
a 



a 



P 



-1 



1 









expl- 



p 



a 



^1 ^ 

v^i2^y j 



^^12^2 



- 1 



P 



-1 



By rearranging terms and simplifying 



F(ai.a2) = 1- 



/ 0 “ ^ 

^12 

. B“ +1 . 
V 12 J 



exp 




v^i2^y j 




noting a-| = B^2^2’ distribution function of the Weibull model 

under bi-axial stress c-\ and 02 can be simplified by 
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a 



F(a,,a2) = 1-exp|-|^-p 

5. From Equation (2-50) 

= J„f,Mf2(v)dudv 









2 1 



2 

noting F-| (oo) = F2(«>) = 1 

^ ' 0 

integrating then we finally got the Equation (2-57) 
R(G-|,a2) = {1-Fi(a)}{1-F2(a)} 
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APPENDIX B 



Structural Applications 



In many structural applications for aircraft, the composite laminae are 
oriented at an angle to the tensile load. Take for example the skin on a wing, 
where the filament angle is oriented primarily along the direction of tension due 
to bending. Transforming the stress by an angle 9 to an off-axis introduces a 
shear coupling that will compensate for the bending/rotational coupling of a 
swept wing. Such a material loading relation is representative of that indicated 
in Figure B-1. We will explore the combined stress state for such a 
configuration and the structural reliability as affected by the off-axis angle 0 . 
Since we are addressing a spatially homogeneous state of stress, both stress 
analysis and the substitution of the combined stress in the probabilistic failure 
criterion can be carried out explicitly. 



X2' 




Figure B-1 : Transformation Of Coordinates 
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For the stress analysis, a load p is applied to the composite material and we 
denote the direction of p to be x-]', then p can be transformed to the principal 
direction of the composite; x-| and X2 where 
xi : fiber direction 
X2 : transverse matrix direction 

and the angle between x-|' and x-|-axis can be defined as 0 which is positive 
clockwise. Using the tensor relation, the stress can be transformed from xi'-x2’ 
axis to x-|-X2 axis as follows: 

cTjj = ami 2nj c^mn' (B-1 ) 

where 

i, j, m, n = 1 , 2 for 2-dimension 

aii = p 

^22’ = ^12’ = ^2l’ = 0 

ajj =^cos(x-|', x-|) cos(xi’, X2)\ =/Cos(0) -sin(0) \ 

^ cos(x2’, x-| ) cos(x2’,X2) ' ' sin(0) cos(0) / (B-2) 

Expanding Equation (B-1) for ajj, 022 > ct 22> then 
<^11 = 3m1 ^n1 crmn' 

= aii an ail' + aii ^21 <^12+^21 ^21 <^22 

= a2i^ ail’ 



an = p cos2(0) 


(B-3) 


For the same reason 




ai2 = -p cos(0) sin(0) 


(B-4) 


G 22 = P sin(0) 


(B-5) 



To calculate the reliability associated with the applied load p, the stress 
components in terms of p (Equations (B-3), (B-4), and (B-5)) can be substituted 
into Equation (2-83) to find the joint failure CDF for the specified external load p. 
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We note that Equation (2-83) is specialized for the case that the effect of 



aj2. snd <522 respective failure modes are independent. Substitution 

yields: 



F(p,e) = 1 - exp 



^pcos^(0)^ 



^psin^(e)^ 
^ ^^2 ) 



a . 



f pcos(e)sin(0) 



I 



P< 



(B-6) 

where 



Oi =aii,G2 = a22,G5 = Gi2 

Comparing Equation (B-6) to (2-83) and (3-18), it can be noted that 
Bi 2 = cot2(0), B-ie = -cot(0) 

Ci2 = Vi2 tan2(0), = -Vie tan(0) 

This joint failure CDF in terms of p and 0 can be plotted in terms of p for an 
specified 0. Similarly, from the Equation (B-6), the linearized failure CDF can be 
obtained by 

psin^(0)| f pcos(0)sin(0) 

(B-7) 

In order to investigate the dependency of the failure mechanism, 
experimental data from Sun and Yamada [Ref. 4] is examined. From the 
experiment, the following parameters were obtained for fiber, transverse, and 
shear force using uni-axial tests, 
tti =20.5, Pi = 127000 
tt2 = 5.66, P2 = 1070 
oc^ = 8.96, Pg = 3078 

Substituting these strength parameters into Equations (B-6) and (B-7), the 




(p,0) = In 



( pcos^(O) 

.l^T 
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Probability Of Failure 



probability of failure for load p applied to different angles can be examined. 
Figures (B-2a) through (B-6b) compare the joint failure CDF to the experimental 
data for different angles. These data will be examined individually for each of 
the five angles where experimental data is available. 




Transverse 
Shear 
Joint CDF 
Experiment 



Figure B-2a : Comparison Of Experimental Data with 
Joint Failure CDF For Independent Case (15 Deg.) 
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Probability Of Failure 



At 0 = 15°, even \when the load is closely aligned with the fiber angle, Figure 
B-2a shows that the fiber strength statistics has no effect on the combined 
failure probability by noting that the fiber failure curve is off-scale for the high p 
region. That is, the failure is always dominated by a combination of matrix 
transfer strength and the matrix shear strength. The interaction of transverse 
and shear can be better observed in the linearized Weibull plot. Figure B-2b, 
from which it can be seen that in the lower tail (low p region) the transverse 
strength dominates and in the upper tail (high p region) the shear strength 
dominates. 




Transverse 
Shear 
Joint CDF 
Experiment 



Figure B-2b : Comparison Of Experimental Data With Joint 
Failure CDF For Independent Case (Linearized, 15 Deg.) 
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Probability Of Failure 



While the experimental data are not close to the predicted joint failure CDF 
under the simplification of independence, its significance can not be evaluated 
because of the small sample size. With the eight samples in the current case, 
neither the upper tail nor the lo\wer tail can be observed,thus the failure 
probability coupling remains unresolved. The location discrepancy perhaps 
has more significance. It suggests that combined transverse and shear 
stresses have weakened the composites and that the failure mechanistic 
coupling needs to be treated. 




Transverse 
Shear 
Joint CDF 
Experiment 



Figure B-3a : Comparison Of Experimental Data With 
Joint Failure CDF For Independent Case (20 Deg.) 
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Probability Of Failure 



At 0 = 20°, Figures B-3a and 3b show that even with a slight angle change of 
5°, the effect of the shear component diminishes rapidly. The transition region 
has shifted higher on the upper tail. Physically, it means that unless a large 
number of samples are tested, only the strongest of the samples will fail in 
shear, and the remaining will fail by transverse stress. For this reason, and 
again because of the small number of samples, the shape indicated by the 
experimental data has no significance. The location difference between the 
data and prediction suggests that mechanistic coupling needs to be included. 




Transverse 
Shear 
Joint CDF 
Experiment 



Figure 3-3b : Comparison Of Experimental Data WithJoint 
Failure CDF For Independent Case (Linearized, 20 Deg) 
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Probability Of Failure 



At 0 = 30°, Figures B-4a and 4b show that the effect of shear is shifted to 
even a higher upper tail. All comments on 0 = 20° data apply to 0 = 30° as well. 




Transverse 
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Joint CDF 
■ Experiment 
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Figure B-4a : Comparison Of Experimental Data With 
Joint Failure CDF For Independent Case (30 Deg.) 
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Probability Of Failure 
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Figure B-4b : Comparison Of Experimental Data With Joint 
Failure CDF For Independent Case (Linearized, 30 Deg.) 
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Probability Of Failure 



At 0 = 45°, Figures B-5a and 5b show that the effect of shear is practically 
nonexistent: in fact, it can not be observed from the range of the scale 
presented. What is significant is the location of data is below the location of the 
joint CDF suggesting mechanistic coupling. 




Transverse 
Joint CDF 
Experiment 



Figure B-5a : Comparison Of Experimental Data With Joint 
Failure CDF For Independent Case (Linearized, 45 Deg.) 
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Probability Of Failure 
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Figure B-5b : Comparison Of Experimental Data With Joint 
Failure CDF For Independent Case (Linearized, 45 Deg.) 
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Probability Of Failure 



At 0 = 60°, Figures B-5a and 5b show that the effect of shear is no longer 
present and that the fit of the data is much improved. The latter observation 
further substantiates the existence of mechanistic coupling. 
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Figure B-6a : Comparison Of Experimental Data With 
Joint Failure CDF For Independent Case (60 Deg.) 
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Probability Of Failure 
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Figure B-6b : Comparison Of Experimental Data With Joint 
Failure CDF For Independent Case (Linearized, 60 Deg.) 
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From the experimental data available in the open literature, it was observed 
that a much larger number of samples are required to identify the probabilistic 
coupling and dependency of strength. What is evident is that mechanistic 
coupling needs to be included as an extension to the investigation herein. 
Furthermore, experimental design using the probability failure criterion is 
mandatory to optimize experimental testing. 
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